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Quantum entanglement, one of the denning features of quantum mechanics, has been demon- 
strated in a variety of nonlinear spin-like systems. Quantum entanglement in linear systems has 
proven significantly more challenging, as the intrinsic energy level degeneracy associated with lin- 
earity makes quantum control more difficult. Here we demonstrate the quantum entanglement of 
photon states in two independent linear microwave resonators, creating iV-photon NOON states as 
a benchmark demonstration. We use a superconducting quantum circuit that includes Josephson 
qubits to control and measure the two resonators, and we completely characterize the entangled 
states with bipartite Wigner tomography. These results demonstrate a significant advance in the 
quantum control of linear resonators in superconducting circuits. 



Quantum superposition and entanglement have been 
demonstrated experimentally using spin-like physical sys- 
tems ranging from atoms to electronic circuits |l|-u|. 
These systems all display strong nonlincarity, and are 
used because this nonlinearity allows straightforward 
quantum control by classical means. The quantum con- 
trol of linear systems, exemplified by the harmonic os- 
cillator, is by contrast more difficult, and has only been 



achieved using nonlinear intermediaries: Atoms 
control optical cavities, ions to control ion motion 
and superconductin g q ubits to control photons in mi- 




crowave resonators [llMl4| . Quantum entanglement of 
cavity photons still presents a significant challenge: Ex- 
periments have demonstrated maximally-entangled pho- 
tons in different polarization modes of the same cavity 
[IH EB| , but the entanglement of photons in two physi- 
cally distinct cavities 17 - 19| has proven more elusive. 

Here we show the deterministic generation of entan- 
gled photon states in two spatially-separated microwave 
resonators, achieved by manipulating the photon states 
with a pair of superconducting phase qubits. We use as a 
benchmark the generation of NOON states [20l - |24l |. com- 
prising a total of N photons in the two resonators (A and 
B), entangled in the quantum state 
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with TV photons in resonator A and zero in B, super- 
posed with the state with the occupation numbers re- 
versed. Such a state has the same degree of entanglement 
as the Bell state, but with N excitations. We also gen- 
erate MOON states, in which, e.g., resonator A has M 
or zero quanta, entangled with resonator B with zero or 
N quanta. We fully characterize the two-resonator photon 
states using bipartite Wigner tomography, which repre- 
sents a non-trivial extension of single-cavity Wigner to- 



(a) 



' — lh^ 


— ii — , 


♦ i 


. 1 

' T 







fj3 



(b) 

6.8 
o>6.4 



6.2 







- \ % - 

. V 

-A— -V 

\ 

ii 


V 

- B -\ 



|0.3 



If}' 

|e)\ 

|g> 




i , — , 6.9 

N 
X 

CD 

o>6.8- 
. 8» 



20 MHz 



0.165 0.185 0.125 0.140 
^bias ( a - u ) 



0.164 0.170 
^bias ( au ) 



FIG. 1: (Color online) (a) Device circuit schematic. Coupling 
resonator C is connected to qo and q\ using 20 MHz coupling 
capacitors, while storage resonator A (B) is connected to go 
(gi) through a 17 MHz coupling capacitor, (b) Qubit spec- 
troscopy, showing probability P e (color bar) vs. microwave 
frequency and flux bias for each qubit. Avoided-level cross- 
ings near 6.8 GHz (dash-dot lines) are due to the coupling res- 
onator C and near 6.3 GHz (dashed lines) due to each qubit's 
storage resonator. Lower right panel shows magnified view of 
circled area, upper right panel shows three qubit levels. 



mography [l|, 0, 12 - 14 1 , and allows us to distinguish en- 
tanglement from an incoherent ensemble. 

To accomplish this goal, we developed a new quantum 
circuit comprising two superconducting phase qubits [25[ 
and three microwave resonators. A sketch of the circuit 
topology is shown in Figs. [Ua). The circuit includes 
a coupling resonator C, connected to both qubits, and 
two state storage resonators A and B, each coupled to 
one qubit. The resonator frequencies arc all different 
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FIG. 2: (Color online) NOON state preparation sequence. 
Resonators are represented by dashed lines, and the qubit 
\g) o |e) and |e) o |/) transitions by dark and light color 
solid lines, respectively. (1): go is excited to \e) and half- 
swapped to C, generating the Bell state |e0) + \g1). (2): 
Coupling resonator swapped to q\, generating \eg) + \ge). 
(3a): N = 1 NOON state ] 10) + |01) generated by fully 
swapping each qubit to its storage resonator. For higher TV 
states: (4): Qubits excited to \ fg) + \gf). (5): One photon 
swapped into storage resonators, generating [eplO) + \ge01). 
Steps (4) and (5) are repeated N — 1 times, generating 
\eg(N - 1)0) + \ge0(N - 1)). (6): Final photon transfer gen- 
erates (7) N = 2 (or higher N) NOON state. 



which allows us to frequency-select the qubit-resonator 
interactions. 

The basic method for generating two-resonator entan- 
gled states, illustrated in Fig. [51 is to excite and then 
entangle the two qubits using the coupling resonator. 
We can swap the resulting Bell state \eg) + \ge) to the 
two storage resonators, creating an TV = 1 NOON state 
1 10) + 1 01) . If we want to generate higher N photon 
states, we instead selectively excite each qubit to its next 
higher energy level |/) (Fig. G^b)), generating the state 
\fg) + \gf), thus using the qubits as "qutrits" 2rj 27|. 
The required microwave excitation is selective, due to 
the anharmonicity of the qubits. The qubit excitation 
is then swapped to each storage resonator through the 
qubit |/) O |e) transition, creating a four-fold entangled 
state |e</10) + |<?e01), where the first two letters indicate 
the qubit states, and the second two numbers the stor- 
age resonator states. We then re-excite the qubits to 
their |/) states, and again swap the excitation to the res- 
onators, generating |eg20) + |ge02). This process can be 
repeated until the entangled state has N—l photons. In 
the final step, each qubit's |e) -O- \g) transition is brought 
on resonance with the corresponding storage resonator, 
swapping the last excitation and leaving the system in 
\ggN0) + \gg0N) = \gg) ® (|JV0) + |0iV», an iV-photon 
NOON state. 

We analyze the final resonator state using the qubits 
as probes. The simplest analysis uses a coincidence mea- 
surement: We bring the qubits into resonance with their 
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FIG. 3: (Color online) Qubit coincidence probability mea- 
surements for iV = 1, 2, and 3 NOON states, and for a mixed 
state. P ge (blue) and P eg (red) oscillate with interaction time 



t at a rate oc y/N, while P ee (brown) always remains small. 
For the mixed state, behavior is nearly identical to the N = 1 
NOON state. Lines are fits to the data. Statistical errors, 
from the measured probability spread of ~ 2 — 3%, are shown 
only for P gg [3(j. Horizontal and vertical axes are same for 
all plots. 



corresponding storage resonators for an interaction time 
r, after which both qubits are measured simultaneously. 
The preparation and measurement sequence is repeated 
~ 10 3 times, yielding the joint-qubit state probabilities, 



Pgg, Pge, Peg, and P ee , where Pge is the probability of 



measuring the first qubit in its ground state with the 
second qubit in its excited state, and so on. We then 
vary the interaction time r, capturing the evolution of 
these probabilities. If a resonator has n photons, the nth 
photon will swap between the qubit and resonator at a 
rate scaling as y/n, while for more complicated states, the 
interaction is a sum of components oscillating at their re- 
spective frequencies, weighted by the photon occupation 
probabilities 13l 14|. 

For resonators entangled in a NOON state, a joint mea- 
surement should correspond to either N photons in one 
resonator and zero in the other, or to the reverse sit- 
uation; the measurement of the qubits "collapses" the 
system onto one or the other outcome. Thus in one 
measurement at most one of the qubits will be in the 
excited state. When averaged over many measurements, 
the maximum probability of measuring a particular qubit 
in |e) is 1/2, while the probability of measuring both 
qubits in |e) should be zero. Therefore we expect that 
P eg and Pge will oscillate between and 1/2, P ee will be 
zero, and P gg should equal 1 — P ge — P eg — Pee- 

Coincidence measurements are shown in Fig. [3] for 
NOON states up to N = 3, and are consistent with 
these expectations. The oscillations seen in these mea- 
surements are however insufficient proof of resonator en- 
tanglement, as an incoherent mixed state can give the 
same results. To demonstrate this, we have controllably 
generated a synthetic ensemble compris ing a 50 % popu- 
lation of 1 10) states and 50 % 1 01) states [30j. Coincidence 
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FIG. 4: (Color online) NOON and mixed state density matrix amplitudes, reconstructed in the photon-number basis from 
bipartite Wigner tomography; states are labeled \mn) where m is the photon number in resonator A and n that in B. Bar 
heights and colors represent matrix element amplitudes. The dominant amplitudes for all three states are in the expected 
locations, although the off-diagonal elements decrease with N, due to the finite qubit dephasing time. Errors for the density 
matrix elements are not shown but are small I30I1. 



measurements (Fig. [3]) of this synthetic mixed state gen- 
erate outcomes identical to those of the N = 1 NOON 
state. 

A more complete resonator measurement, that can re- 
solve entangled from mixed states, uses bipartite Wigner 
tomography, a nontrivial extension of single resonator to- 
mography 9 



13l . Il4j . This involves injecting a coherent 



Gaussian microwave pulse into each of the storage res- 
onators, with controlled amplitude and phase, displacing 
the resonator states in phase space. The resonators are 
then simultaneously measured with a joint probability 
measurement, now as a function of the amplitude and 
phase of the coherent pulses. From the complete set of 
measurements, the two-resonator density matrix can be 
calculated [3(ij |. 

In Fig. S] we display the amplitudes of the density 
matrices measured for resonator NOON states up to 
TV = 3, as well as for the mixed state. While there 
are non-idealities, the desired non-zero matrix elements 
are clearly apparent for the NOON states, while for the 
mixed state, the density matrix has only zero- valued off- 
diagonal elements. The state preparation fidelities, F = 
(i/j\p\ip), are found to be 0.76 ±0.02 (N = 1), 0.50 ±0.02 
(N = 2), and 0.33 ± 0.02 (N = 3). For N = 1, the most 
probable entanglement of formation [28| is EOF = 0.51, 
while for N = 2 and N = 3, EOF = 0.31 and 0.28, re- 
spectively; for the mixed state, the EOF is zero. We also 
calculate the negativity N e (p) — J2j max (0, Hj), where 
Hj are the eigenvalues of the partial transpose p PT of the 
density matrix 22, 29|, and N e (p) > indicates entan- 
glement. The negativities are found to be 0.56 ± 0.03 
(N = 1), 0.32 ± 0.03 (N = 2), and 0.27 ± 0.01 (N = 3); 



for the mixed state we find zero with an upper bound of 
0.001. The decrease of these values with photon number 
N is compatible with expectations: The state prepara- 
tion requires phase-coherence of the four-element entan- 
gled states for most of the preparation sequence, which is 
limited by the qubit coherence time T 2 |30|. Other than 
this technical limitation, the deterministic generation is 
completely scalable to large N. 

A hallmark of NOON states is their rapid phase evo- 
lution [2^, [21], 23, 24 1, which can be verified by Wigner 
tomography using two distinct methods. For the N = 1 
state, after entangling the qubits in a Bell state, we wait 
for a variable time and then swap the state into the stor- 
age resonators. The density matrices measured at three 
different delay times are shown in Fig. [5ja). The phases 
of the off-diagonal elements rotate with time due to the 
qubit-resonator frequency difference, as in Fig. EJb), 
showing the expected linear dependence. 

This phase-measurement method suffers from the short 
qubit dephasing time. A second method is to change 
the phase reference for the coherent pulses used in the 
Wigner tomography, avoiding storage of the state in the 
qubit. We add an additional phase to the pulses applied 
to resonator A only. The resulting density matrices show 
the expected rotation of the off-diagonal elements. In 
Fig. [5jc) we plot the off-diagonal phase angle for different 
N; the N = 3 state evolves three times faster than the 
N = 1 state, as expected. 

We also used tomography to measure the NOON state 
decay [30]. We find that the off-diagonal elements de- 
crease at approximately the same rate as the diagonal 
elements, with a decay time tjj « 3 /is, consistent with a 
Markovian environment 14 1 . 
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FIG. 5: (Color online) Phase sensitivity of NOON states up 
to N = 3. (a) Density matrix for N = 1 NOON state at dif- 
ferent times after qubit entanglement. Each element is repre- 
sented by an arrow, with orientation determined by the phase 
angle (scale on bottom left). Off-diagonal elements rotate 
with delay time (black: 0ns, blue: 12ns, red: 24ns), due to 
frequency difference between the qubit operating point and 
resonator. Dephasing causes decrease in amplitude of off- 
diagonal elements, (b) Phase angle of the upper left non-zero 
off-diagonal element in (a) versus delay time. Line is fit to 
a rotation rate of A/ = 21.6 ± 0.8 MHz, which corresponds 
well to set frequency difference. Error bars indicate maxi- 
mum phase uncertainty, (c) Rotation angle of off-diagonal 
element versus controlled, additional phase angle used in co- 
herent state displacement. Lines are fits (slopes indicated on 
plot), consistent with expected phase sensitivity. Confidence 
bounds are as in (b), with uncertainty increasing with N due 
to increased phase sensitivity. 



We note that the generation sequence allows an ad- 
ditional flexibility: We can add different numbers of 
quanta to the resonators, thus generating MOON states, 
|t/>) = I MO) + IO/V). An example with M = 2 and N = 1 
is detailed in [30|. 

The capabilities we have demonstrated here, generat- 
ing complex entangled photon states in two resonators, 
hold promise for new quantum architectures in which 
superconducting resonators play a more central role. 
The protocol can be extended in a natural way to 
entangle larger numbers of resonators, allowing, e.g., 
the direct generation of resonator GHZ and W states 
@, S B EH, HH HH • The longer coherence times achievable 
in superconducting resonators will be of direct utility in 
performing more complex quantum algorithms, further- 
ing the capabilities of superconducting quantum circuit 
architectures. 
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Supplementary Material for 
"Deterministic entanglement of photons in two 
superconducting microwave resonators" 

MATERIALS AND METHODS 

The device fabrication is similar to that published pre- 
viously [s"H . The half- wavelength superconducting copla- 
nar waveguide resonators are made of rhenium deposited 
on a c-axis single-crystal sapphire substrate, with a 5 /zm- 
wide center signal trace and 10 ^m gaps to the ground 
plane metallization on either side of the center trace. We 
place a single lithographed shorting strap connecting the 
two ground planes at the midpoint of each resonator to 
improve the quality of the grounding. This point is a volt- 
age node for the fundamental half-wave resonant mode, 
so that there is minimal additional dielectric loss from 
the shorting strap's underlying amorphous Si insulating 
film [S3. 

In the circuit layout, the coupling resonator C is de- 
signed to have a higher resonance frequency than the 
two state storage resonators A and B. This prevents the 
qubit frequencies from having to cross the C resonator 
frequency during NOON state amplification. The two 
storage resonators A and B are designed with slightly 
different resonance frequencies, to avoid possible inter- 
ference between the resonators. The full frequency span 
in the design was chosen to be about 550 MHz, within 
the dynamic range of our custom microwave electron- 
ics. The two superconducting phase qubits and copla- 
nar waveguide resonators are fabricated together, using 
our standard multi- layer process S3]. We use interdigi- 
tated coupling capacitors between the qubits and the res- 
onators, calculated to each have a coupling capacitance 
of 1.9 fF. The actual coupling strengths vary slightly with 
resonator frequency; the detailed component parameters 
are listed in Table [Si] 
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FIG. SI: (Color online) One-photon swaps between qubit go 
and resonator A. (a) Pulse sequence for measuring swaps be- 
tween the qubit |e) <H> /) transition and resonator A. When 
off-resonance, the qubit is excited by a pair of n pulses, with 
the first pulse tti taking the qubit from \g) to |e) and the sec- 
ond pulse 7T2 making the transition from |e) to |/}. The qubit 
e) <-> I /) transition is then tuned close to resonance with res- 
onator A and the qubit and resonator left to interact for a 
time r. The triangular pulse at the end is used to measure 
the probability of the qubit being in the |/} state at the end of 
the sequence, (b) Left: Qubit \g) <-» |e) swaps with resonator 
A. The qubit is prepared as in (a), but without the second 
pulse 7T2. The qubit |e) probability (color bar) is plotted as a 
function of interaction time r and frequency tuning $bias- Os- 
cillations in probability are due to the qubit excitation swap- 
ping with the resonator. Right: Qubit |e) -o |/) swaps with 
resonator A; the preparation of the qubit is as in (a) and 
the plot shows the |/)-state probability as a function of in- 
teraction time r and qubit flux bias ^bias- Red dashed lines 
indicate the on-resonance "^bias values used in the experiment. 
The swap frequency for the on-resonance je) \f) transition 
is 1.403 « V2 times that for the on-resonance \g) «-> \e) tran- 
sition. 



GENERATION SEQUENCE TUNE-UP 

The time required for each qubit-resonator i-SWAP 
is calibrated separately. The swap times obtained from 
these calibrations scale correctly as * s fn with the number 
of photons n in the resonator [S5j and also depend on the 
state of the qubit. Examples of the swap calibrations for 
a one-photon swap between qubit qa and resonator A are 
shown in Fig. EU for both the \g) O |e) and |e) O |/) 
transitions. The swap time for the |e) 4-> \f) transition is 
approximately 1/V2 times that for the \g) -H> |e) transi- 
tion. This scaling is as expected, as the multi-level phase 
qubit can be well-approximated as a weakly nonlinear 
harmonic oscillator for the energy levels confined by the 
qubit's metastable potential well. The scaling confirms 
that we can use harmonic-oscillator-like raising and low- 



ering operators for the three-level qubit, as was done in 
the numerical simulations (see below). 

Figure [S2] shows the detailed pulse sequence used to 
generate and measure the N = 2 NOON state. Sequence 
steps in general are calibrated and checked separately to 
maximize preparation fidelity, when possible. For exam- 
ple, we first optimize the qubit Bell state preparation at 
the end of step I in the preparation sequence detailed in 
Fig. [S2l The fidelity for the Bell state is above 0.80, with 
entanglement of formation 0.59, which agrees well with 
numerical simulations, performed using a pure dephasing 
time T<f, = 300 ns for the qubits (see below) . State tomog- 
raphy of the qubits is also as expected: Fig. IS27 b) shows 
the density matrices extracted from coupled qubit to- 
mography, measured at different times during the N = 2 
NOON state preparation. We note that at the end of 
the sequence, both qubits should return to their ground 
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TABLE SI: Resonator and qubit parameters. The qubit parameters are quoted for when the qubits are off-resonance (see 
below). The qubit nonlinearity /nonlinear is the frequency difference between the \g) o |e) and \e) -H- \f) transitions. The 
phase coherence time T^> is obtained using a Ramsey interference experiment, which yields the Ramsey time Ti, from which 
we calculate l/T^ = I/T2 — l/(2Ti). measured for resonators similar to those used here [s3 ] satisfies T$ 3> Ti; we assume 
the same applies here. T$ for the qubits decreases with increasing length of the pulse sequence due to the 1// nature of the 
phase noise. For most state generation sequences used in this experiment, of typical length 50 to 100 ns, the qubit is in the 
range of 200 to 300 ns. Coupling strengths correspond to the measured splitting in frequency units, with g appearing in the 
Hamiltonian shown in Eq. ()S1|I . 



states. Experimentally we observe small populations in 
the excited states |e) due to decoherence and pulse im- 
perfections. The exact qubit state after the NOON state 
generation is measured and is used as the initial state for 
the qubits when performing Wigner tomography on the 
storage resonators (see below). 



NUMERICAL SIMULATIONS 

Numerical simulations were performed using the model 
Hamiltonian 

H = E H « + E ft*' [ b t h i + \ 
+ E n 9n<*3 ( a i b j + a « 6 i~) 

+ E ^ ["«(*)<+"« (*)*««]> ( S1 ) 



where for simplicity we assumed a constant nonlinear- 
ity f\g)++\e) - f\e)**\f) = 200 MHz, so that f\g)^\f) ^ 
2/| g )<->.|e) — 200 MHz. We approximated the multi-level 
qubit a qi and a qi by the raising and lower operators for 
the lowest three levels of a harmonic oscillator, as dis- 
cussed above. 

Decoherence was approximated using the Lindblad 
master e quat ion taking into account the Markovian envi- 
ronment [S6|, where two characteristic decay times, the 
energy relaxation time T\ and the pure dephasing time 
Tfj), were used for each resonator and qubit. 

The simulations do not directly account for the non- 
Markovian 1// character of the phase noise in the qubits. 
To account for this, we used a sequence-time dependent 
for each qubit, as obtained from Ramsey interference 
measurements. We used T$ = 300 ns for ~50 ns-long 
sequences and T$ = 200 ns for ^100 ns-long sequences. 
The resulting simulations agree reasonably well with the 
experimental measurements. 



where H qi is the Hamiltonian of the qubit qi, a+ and 
a qi (b~j and bj) are the raising and lowering operators 
for the 3- level qubit qi (resonator j), g qi ^j is the cou- 
pling strength between qubit qi and resonator j, with a 
sum over all possible qubit-resonator combinations, such 
that ( qi ,j) G {(q ,A),(q ,C),(q 1 ,B),(q 1 ,C)} and Cl qi (t) 
is the time-dependent, two-tone (/| ff wi e \ and /i e )++|/>) 
microwave drive on qubit q^. 

The 3-level qubit Hamiltonian H qi was approximated 

as 



H„ 





hf\ g )++\ e ) 
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BIPARTITE WIGNER TOMOGRAPHY 
Displacement Pulses 

The bipartite Wigner tomography is an extension of 
a method described elsewhere js?], [sB], S4|. We dis- 
place each resonator with Gaussian pulses \a) and \j3) 
(resonator A and B, respectively), with variable phase 
and amplitude. The values of a and j3 are distributed 
over several concentric circles in the complex plane, cen- 
tered on the origin, where the distribution of values 
varies approximately with the size N of the NOON 
state. The radii of the circles run through the set 
rj G |0, 0.2, 0.7, 0.9, 1.3}, in square-root of photon num- 



ber units 



The pulse values are evenly distributed on 
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FIG. S2: (Color online) NOON state preparation, (a) Sequence to generate and measure TV = 2 NOON state. Qubit transition 
frequencies \g) -H> |e) (dark-color solid lines) and \e) <~> j/) (light solid lines) are tuned to control coupling of each qubit to 
resonators (dashed lines); qo in red, qi in blue. Vertical lines divide sequence into steps; circled numbers indicate points where 
qubit tomography was performed. Step I: A 7T] g )_j.| e ) pulse excites q to |e), followed by a V i-SWAP entangling qo and resonator 
C. A full i-SWAP between C and q\ disentangles the resonator and generates the qubit Bell state \eg) + \ge). For the TV = 1 
NOON state this entanglement is transferred to the storage resonators (not shown). For TV > 1, Tri,,)^^ pulses excite each 
qubit to yield \fg) + \gf). Both qubits' je) -H> /} transitions are then tuned to the respective storage resonators, and an 
i-SWAP transfers one photon to the resonators. For TV = 3, this is repeated to transfer a second photon (not shown). In 
step III, both qubits' |<j) o |e) transitions are tuned to the resonators and an i-SWAP transfers the remaining excitation to 
the resonators, completing the generation, leaving both qubits in their ground states. Step IV: State analysis. The storage 
resonator A (73) is displaced in phase space by a Gaussian pulse \a) (|/3)), followed by an on-resonance interaction with each 
qubit's \g) o |e) transition for a time r, followed by joint-qubit state readout. Many repetitions of this sequence are combined 
to generate joint-qubit state probabilities P gg , P ge , P eg , and P ee versus r. The joint-qubit state probabilities can be used to 
obtain the two-resonator density matrix (see text), (b) Qubit density matrices extracted from tomography at circled points 
in (a). (1) Qubit Bell state. (2) Qubit state after transferring first photon into resonators. This is a state with qubits and 
resonators entangled, so the density matrix is as expected. (3) Qubit state after TV = 2 NOON state generation. Both qubits 
are very nearly in their ground states, and thus can be used to measure the resonators. 



each circle, with complex values Tj exp(i2ir£/N rj ), where 
£ ranges from 1 to TV r . and TV r . is an integer ranging from 
1 (rj = 0), 5 or 6 (rj = 0.2), up to 9 to 15 (r,- = 1.3), 
depending on the number of photons in the NOON state. 

We use every possible combination of values of a and (3 
distributed over the circles of the same radius for tomog- 
raphy, i.e., for each value of a, we use all (3 values with the 



same amplitude as a. The total number of displacement 
pulse combinations is thus quite large and increases with 
the photon number TV in the NOON state, typically in- 
volving of order a few hundred pulses. The displacement 
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pulses can be expressed as 
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By bringing both qubits (initially in their \g) states) on 
resonance with the resonators, the joint number states 
contained in two resonators, i.e., the diagonal elements 
of p, can be read out, as each diagonal element swaps 
with the qubits at a different rate [S5|, resulting in a 
distinct time-dependence for the probabilities P gg , P ge , 
P eg , and P ee . These can be numerically simulated using 
the device parameters from Table IS11 

As displayed in Fig. IS2f b). there is a small non- 
zero occupation of the excited state of each qubit af- 
ter the state generation sequence, due to decoherence 
and pulse imperfections. We use the measured qubit 
state after the state generation sequence as the qubit 
initial condition when numerically simulating the tomog- 
raphy. Using these simulations, we obtain the time 
dependence for the each of the probabilities P gg , P ge , 
P eg , and P ee corresponding to different combinations 
of photon number (Fock) states in the two storage res- 
onators. Examples of these probability traces are shown 
in Fig. IS3f a) for some selected initial states. The time- 
dependent traces for these probabilities, for the set of 



Fock states {|m) 



,m = 0, l,2,...,n = 0,1,2,...}, 



are then used to decompose the experimentally-measured 
time traces, which yields the probability distribution for 
the Fock number states contained in the storage res- 
onators. This thus yields the diagonal elements of the 
experimentally-measured displaced density matrix p. 

We obtain the diagonal elements of p by doing a least- 
squares fit of the time-dependent probabilities, corrected 
for measurement fidelity. We use the MATLAB pack- 
ages YALMIP and SeDuMi for the fitting. The num- 
ber of fitting parameters is the number of diagonal el- 
ements, directly determined by the maximum photon 
number state contained in the resonators, plus the num- 
ber of photon quanta added by the displacement pulses 
\a) and \/3). Fits are done with constraints P mn > and 
J2 m n Pmn = 1 to return meaningful probability values. 



Examples of these fits are shown in Fig. 
NOON state. 



for the N = 1 



TABLE S2: Density matrix elements as a function of delay 
time for the N — 1 NOON state shown in Fig. [S4] Uncer- 
tainties for the diagonal elements are in parentheses, while 
the magnitude of the uncertainty for the off-diagonal terms 
are given in the last column. 



Joint-Resonator Density Matrices 

With the diagonal elements of p measured for a set 
of displacements {Dab(&j ■, Pk)}, /"initial can be solved 
for by inverting Eq. (|S4[) through a linear least-squares 
fit, while restricting the density matrix to be Hermitian. 
The resulting density matrices may have small negative 
eigenvalues due to noise. We use the MATLAB packages 
YALMIP and SeDuMi to perform semi-definite convex 
optimization programming, allowing us to find a physi- 
cal density matrix that is closest to the actual matrix. 

When solving for pinitiai of a NOON state, we restrict 
the dimension of pinitiai to an N x N matrix, even though 
the dimension of p can be significantly larger than this 
due to the displacement pulses. We zero-pad the ele- 
ments in pinitiai that have photon indices larger than N. 
This approach is validated by the coincidence measure- 
ments (see main text), as we do not detect any frequency 
components for number states above N prior to injecting 
a displacement pulse. 



NOON STATE DECAY DYNAMICS 

The Wigner tomography allows us to study the deco- 
herence dynamics of the bipartite system [SJ]. The ex- 
perimental results, compared with numerical simulations, 
are shown in Fig. IS41 with relevant elements in Table [S2| 
We note that the time evolution of the off-diagonal ele- 
ments in the two-resonator density matrix, which repre- 
sent inter-resonator coherence, is different from the evo- 
lution of the corresponding off-diagonal elements for a 
single resonator, which represent intra-resonator coher- 
ence [S4|. 
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FIG. S3: (Color online) Joint-state probabilities and photon number distribution analysis for N = 1 NOON state, (a) Joint- 
qubit state probabilities from numerical simulations, for a few selected initial resonator states in A and B, as indicated. Lines 
are a guide to the eye. Device parameters are as in Table [Si] (b) Left panel, measured joint-qubit state probabilities for null 
displacement pulses a, ft = 0. Lines are fits to the data. The fit joint-photon probability distributions P\m) A \n) B are displayed 
in the right panel. The most significant occupations are in |0}a|1)s and |1)a|0)b, as expected for this state, (c) Left panel, 
measured joint-qubit state probabilities for the displaced state with a — 0.7 and /3 = 0.7 (square- root of photon number units). 
Lines are fits to the data. The fit joint-photon probability distributions P\m) A \n) B , shown in the right panel, are significantly 
different from those shown in (b) due to the displacement. The matrix elements for the higher photon number states are 
omitted for viewing, although these may not be negligible for larger displacement pulses. 



Error Analysis 



Validation 



Statistical errors in the qubit probability measure- 
ments as well as uncertainty in the amplitude calibra- 
tion for the displacement pulses ctj and (3u are used to 
estimate the uncertainty in the amplitude and phase of 
each density matrix element. These errors are found to 
be small, in part because the constraints on the analysis 
filters unrealistic values. Instead we find that slow phase 
drifts in the electronics, perhaps dominated by ambient 
temperature fluctuations, give the main phase uncertain- 
ties, especially during long measurements. Evaluating a 
single density matrix usually takes a relatively short time 
during which these drifts are minimal. However, measur- 
ing a series of density matrices such as Fig. 5 in the main 
paper, takes a much longer time, typically 12 to 20 hours, 
allowing for more significant drifts. These mostly affect 
the phases of the density matrix elements, rather than 
the amplitudes. 



The bipartite Wigner tomography was validated by 
several consistency checks. (1) The density matrix is 
as expected for a range of different states, including 
the highly entangled NOON states, the energy eigen- 
states |2)ji|0)b (Fig. !55Ta)). the separable (product) 
state (|0) - |1»a(|0) - |l))s (Fig. ISgTb)) and the un- 
cntangled mixed states (see below). (2) The NOON 
states display the expected phase sensitivity as a func- 
tion of photon number N, as shown in Fig. 5 in the main 
paper and Table [HI (3) The NOON state fidelity and 
entanglement of formation agree reasonably well with nu- 
merical simulations. (4) The time-dependence of the den- 
sity matrix elements, showing uniform exponential decay 
of all elements, is as expected and agrees with numerical 
simulations, as shown in Fig. IS4I and Table [S2l (5) The 
calculated negativities are significantly above zero for the 
NOON states and precisely zero (within the measurement 
error) for the unentangled mixed state (see below). 

We note that bipartite Wigner tomography can mea- 
sure any matrix element with a relatively high accuracy, 
as we can displace the system by an arbitrary amount 
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FIG. S4: (Color online) Decay of the N = 1 NOON state, (a) Density matrix amplitudes for the N = 1 NOON state at 
selected delay times. Bar heights and colors represent the matrix element amplitude, scale on right. The off-diagonal element 
amplitudes decrease at a rate close to that of the diagonal counterparts. Values for the relevant elements are tabulated in 
Table [S2l (b) Left: Experimental data, showing amplitude of the off-diagonal element as marked by a green ellipse in the first 
panel in (a), plotted versus post-preparation time. Line is a fit yielding a single decay time td, consistent with a Markovian 
dissipative environment [S8j |. Right: Numerical simulation of the same density matrix element as the experiment, with the line 
an exponential fit. The simulation uses the measured density matrix as the initial condition and the resonator Ti value listed 
in Table [ST] The simulated decay time to is in good agreement with experiment. 
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FIG. S5: (Color online) Benchmark resonator states used to validate the bipartite Wigner tomography, (a) Photon-number 
representation of the two-resonator density matrix (error bars not shown), for the energy eigenstate |2)a|0)s, generated by 
pumping two photons into resonator A while leaving B in its ground state. The magnitude of each matrix element is represented 
by the full length of the corresponding arrow and the phase angle determined by the direction of the arrow in the complex plane 
(scale on bottom left). Density matrix is as expected from state preparation, (b) Two-resonator density matrix (error bars not 
shown) for the test resonator state (|0) — |1))a(|0) — |l))s, generated by creating the \g) — i\e) state in each qubit, followed by 
a complete i-SWAP transfer to each storage resonator. Non-idealities are likely due to inaccuracies in the preparation pulses. 



in phase space. Even for relatively small off-diagonal ele- 
ments, such as the desired off-diagonal term in the N = 3 
NOON state, the tomography can unambiguously evalu- 
ate this element and measure its sensitivity to external 
phase perturbations. 



Ensemble of Mixed States 

We use a synthetic ensemble of mixed states to illus- 
trate the hazards involved in relying purely on coinci- 
dence measurements for demonstrating NOON-state en- 
tanglement. The ensemble comprises a mixture of 50% 
|1)a|0)b and 50% lO)^!!)^ states, i.e., an ensemble with 
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TABLE S3: (Color online) Density matrix elements for the 
data shown in Figs. 5(a) of the main paper, coded by color 
(uncertainties not shown). The phase uncertainty in the off- 
diagonal elements is dominated by the slow phase drift of our 
electronics (see text). 

the same probability of being measured in either of the 
states forming the N = 1 NOON state, but without any 
entanglement. This is done by generating the pure state 
|1)a|0)_b and measuring the time-dependent joint proba- 
bilities P gg (r), P 9 e(r), P eg (r), and Pee(i~) for this state. 
We then generate the other component of the ensemble, 
the pure state |0).a|1)b, and repeat the generation and 
measurement procedure. Each value of r involves 300 re- 
peats of the preparation and measurement sequence for 
each of the pure states. We then combine the measure- 
ment results with equal weights, creating the joint prob- 
abilities for the synthetic ensemble; these data are shown 
in the main paper. The tomographic analysis yielding 
the density matrices is done in the usual way. The out- 
come of the ensemble measurements are shown in the 
main text, with the joint probabilities evolving in a way 
indistinguishable from the entangled NOON states, but 
the density matrix for the ensemble revealing a complete 
lack of entanglement, as witnessed by the negligible val- 
ues for the off-diagonal elements. 



A, yielding the state \eg(M - 1)0) + \ge0(N - 1)). The 
final qubit excitations are then transferred in step III, 
resulting in the MOON state \gg) ® (|M0) + |0iV>), with 
the qubits disentangled from the resonators. A MOON 
state generated in this fashion, with M = 2 and N = 1, 
is shown in Fig. IS6l 
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FIG. S6: (Color online) Measured density matrix (errors not 
shown) of the MOON state with M = 2 and N = 1. The 
generation sequence is as described in the text. The state 
fidelity is 0.42 ± 0.01, lower than our typical NOON state 
fidelity, due to technical issues with the particular device used 
for this experiment; the entanglement of formation is 0.16. 
The negativity is N e — 0.14 ± 0.01, indicating a statistically 
significant entanglement. 



MOON STATE 

The NOON-state generation protocol can be simply 
generalized to generate MOON states, with different pho- 
ton numbers in the two entangled resonators. The gen- 
eration is similar to the NOON state sequence shown in 
Fig. IS2I We assume M > N: After generating the Bell 
entanglement between two qubits at the end of step I 
(Fig. IS2fa)). we repeat step II N — 1 times, yielding the 
four-fold entangled state \eg(N- 1)0) + \ge0(N- 1)). The 
photon amplification and transfer process (step II) is then 
applied M — N times, but only to qubit qo and resonator 
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